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Mixed-Path Ground- Wave Propagation: 
2. Larger Distances 

James R. Wait and James Householder 

The theoretical results given in part 1 (NBS Research Paper 2687) for ground wave 
propagation over a mixed path on a flat earth are generalized to a spherical earth. The 
problem is formulated in terms of the mutual impedance between two vertical dipoles which 
are located on either side of the boundary of separation. Extensive numerical results are 
given in graphical form for a mixed land-sea path at frequencies of 10, 20, 50, 100, and 200 
kilocycles per second. 

1. Introduction 

r 

In a previous paper [1] l (designated hereafter as part 1), calculations were presented for 
the propagation of vertically polarized ground waves across a boundary separating two con- 
trasting homogeneous media. The earth was assumed to be flat so that the results were re- 
stricted to short distances. It is the purpose of this sequel to extend the analysis and computa- 
tions to larger distances, where the curvature of the earth must be considered. 

The mutual impedance, Z', between two vertical electric dipoles located at A and B, de- 
picted in figure 1, is considered. The earth medium to the left of the boundary line has a 




Figure 1. Model for -propagation between A and B straddling a land-sea boundary on a smooth spherical earth. 



conductivity, a, and dielectric constant, e, and the corresponding constants for the medium to 
the right of the boundary line are a x and e^ Using the same approach as in part I, it follows 
that 



Z ' =Z+ A!I ["\xa;-3xhJ *. 



(i) 



where Z is the mutual impedance between the dipoles over a homogeneous curved surface of 
conductivity a and dielectric constant e. E A , E&, H A , H' B are the electric and magnetic fields 
tangential to the surface S which is the portion of the curved surface to the right of the boundary 
line. The unprimed quantities refer to the case where the dipoles are in the presence of a 
corresponding homogeneous surface of constants a and e throughout. The primed quantities 
refer to the case where A and B are on the modified surface. The subscripts A or B on E A , E' B , 
etc., specify that the source of the respective field is a dipole carrying a current, /, at A or B, 
respectively. 



Figures in brackets indicate the literature references at the end of this paper. 
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Invoking the principle of stationary phase, the surface integral in eq (1) can be reduced to 
a line integral between the coast line at O and the dipole B. The reasoning here is almost 
identical to that used in part 1. AOB, rather than being a straight line on a flat earth, is now 
the path of stationary phase between A and B. In the case of a spherical earth, AOB is on a 
great circle. 2 In the present problem, the concept of surface impedance is introduced as it was 
in part 1, the tangential electric and magnetic fields are assumed to be related by a complex 
constant of proportionality. For the medium to the left of the boundary, this constant is 
t] and the corresponding constant for the medium to the right of the boundary is rj u where, in 
ohms, 



and 



„=120^(l+g) H (2a) 

»=120^(hD* (2b) 



with /3=27r/wavelength, a=(iafji,o>)— eixu 2 ) 1 / 2 , a 1 =(ia 1 fjLco—eifxo)Y /2 and /*=47rX10~ 7 h/m. 

The mutual impedance, Z, between the short antennas A and B for the homogeneous ground 
of electrical constants a and e is given by 



hghbiya) 
2ird l 



-^W,,,a)[l+^^ 2 } (3) 



where d is the great circle distance AOB, h a and h b are the effective heights of the dipole anten- 
nas at A and B, and W(d,7),a) is a slowly varying function that approaches unity for a flat, per- 
fectly conducting earth. The function W is obtainable from the Van der Pol-Bremmer theory 
[2] for the field of a dipole over a homogeneous sphere. It is given, in convenient form, by 

«'<«.<.«K?)" V "' 4 <«"' § "^y (4a) 

with 5=-i(0a)- i/3 (l2OT/r 1 ) and (4b) 

1 1 + 4^ 



2v,8 SvlS 2 12^5 3 



(4c) 



and 



v =0.S0Se- iir/3 
v l =2.577e~ iT/z 

j, 2 = 3.824e-* r/3 

^^l/2[37r(s + l/4)PV i7r/3 for s>2. 

An expansion for the coefficients t s in a power series in 5 is also available [2]. In the preceding 
formula, a is the effective radius of curvature of the surface in the principal plane containing 
AOB. 

The mutual impedance for the inhomogeneous path is now written in the form 

where W is some function of d, rj, rji, and a and can be expected to be slowly varying compared 
to e~ m and must necessarily approach unity as d tends to zero or if rj and 771 tend to zero and a 
approaches infinity. 



2 The coastal refraction effect is neglected here because for most cases it is very small, as shown in part 1. 
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The essential modification here to the flat-earth theory in part 1 is to replace the functions 
F(d,r]) and F'(d,ri,iii) by W'(Y//r?,a) and W'idtf^a). It then Follows that 



"■^wW^^S)!* NHf** ^ 



(6) 



which is the spherical-earth counterpart of eq (13) of part 1. In the above, a, the integration 
variable, ranges from B to O (i. e., rfi = OB). This result is an approximation that is not valid 
near the boundary line because it utilizes a stationary phase principle, as in part 1. Further- 
more, the wave reflected from the boundary is neglected. Actually, this effect is very small at 
distances greater than a wavelength from the boundary line. 

The prime task now is the evaluation of the integral in eq (6). Inserting the Van der 
Pol-Bremmer formula for the W's, the integrations can be carried out to yield 

/ _ x - - exp[~ir s (fia) u:i dla}[exp[i(r-r])^ (*a)»«"|-l"|, 

(7) 

where the coefficients r s are given the eq (4c) et seq. and t\ are identical in form to r,, except 
that o-j and e { replace a and e. This double summation is very cumbersome and converges 
poorly. An alternative formula may be derived by employing the following series formula 
[3, 4] for the W's: 

/rf — a \n/2 

where 1/T= — (?'/3/2)(r7/1 20tt) 2 . The, coefficients up to order 9 are given by 
A =l, i4,= -iV5r, A 2 =-2, A 3 =i^(l-5 3 /2) 

A,=\ (1-5 3 ), A=-^(l-35 3 /2), A 6 =- ] | ? (l-25 !( +|a 6 ), 

where 5=—i(pa) 1/ * (\2{)w/r)). The corresponding expression for W{a y 7] l} a), for example, is iden- 
tical to the above with d— a y T, and 5 being replaced by «, T x and 8 l} respectively, where 

1/Ti= -(i/3/2)(V!207r) 2 and 5i = -i(0a) 3 (120^). Therefore, 

m",ma)=£Al(£) n ", (10) 

m=0 \J-1/ 

where the A} n coefficients are functions of 8 { . After carrying out the integration indicated in 
eq (6), the formula for W takes the form 

WUv,m,a)=W(d, v ,a)-(pJX r ^P) S £ A n AU m , „ (11) 

W here 

r<h m zl n ~ 1 

/»,„=J o « 2 ((*-«) 2 ^. (12) 
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(9) 



The integral I mt n can be reduced by the formulas 



m + lK 2 ' " 3 > 



and 



/ -= 2tan " 1 (A) M 



The preceding expression for W r/ is useful when distance d is small compared to T and TV 
In view of the fact that extensive numerical data for the W functions are available [5, 6] r 
it would seem to be preferable to evaluate the integral by a numerical method. Certainly this 
appears simpler than attempting to sum. the doubly-infinite series formulas discussed above. 
To convert the integral to a form, suitable for numerical integration, it is desirable to change 
the variable of integration to remove the singularity at a=0. Then, letting a=x 2 , it follows 
that 



W '=W+2 M A \ l -^\ e-n P' W(d-Av,a^Y^' a)dx , (14) 



where 



a + icoe 2 V1207T/ 



At low radiofrequencies, displacement currents in the ground are small, and consequently k 
and po are real and positive. 

For purposes of computation, the following values of the parameters are employed: 

a=lQ~ 2 mho/m corresponding to land, 
0^ = 4 mhos/m corresponding to sea, and 
a = 4/3 times the earth's radius. 

In figure 2 the amplitude of W is shown plotted as a function of distance d for frequencies of 
10, 20, 50, 100, and 200 kc. Two values of d h ld are indicated on the curves, namely, cor- 
responding to an all land path, and 0.8, corresponding to an 80 percent sea path. For the lower 
frequencies, the curves are almost indistinguishable, and the shape of W is determined pri- 
marily by diffraction by the earth's curvature. The phase lag <£'(i. e., $'= — arg W) is shown 
plotted in figures 3 to 7 as a function of d for the same frequencies. The values of djd, which 
indicate the ratio of the length of the sea path to the length of the total path, are 0, 0.2, 0.4, 0.6, 
0.8, and 1.0. 

The preceding analysis has been discussed in relation to the mutual impedance between two 
dipole antennas at A or B. It is probably desirable to express the results in terms of the vertical 
electric field, E, at B (or A) for a standard source at A (or B). For example, 

E=EoW'=EoW\e- i +', (15) 

where 



d L my Pdj ' ' 



o e -W. (16) 
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Figure 2. Amplitude factor for mixed-path 
propagation. 



Figure 4. Phase factor for mixed-path propagation. 
f=20 kc. 
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Figure 3. Phase factor for mixed-path propagation 

/=10kc. 
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Figure 5. Phase factor for mixed-path propagation. 
/=50kc. 
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Figure 6. Phase factor for mixed-path propagation. 
/=l00kc. 
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Figure 7. Phase factor for mixed-path propagation. 
/=200 ke. 
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Figure 8. Amplitude for near field. 




I'i < i re 10. Phase of the ground wave for propagation from land to sea, 
(Transmitter Is distance d> from the coast) /=100 kc. 
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Figure_9. Phase correction for near field. 
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The "near field", E , is chosen such that the radiation component is 0.10 v/m at d= 1 mile= 1,605 
m. The amplitude, E , and the phase, <£ , are shown plotted in figures 8 and 9, respectively, for 
frequencies of 10, 20, 50, 100, and 200 kc. It should be noted that |2? | = 160.5/d and $ ^ 1/pd 
when j3d^>y>l. Actually, E is the total field of the source on a flat, perfectly conducting earth. 
To a good first approximation the total field E over a mixed path on a spherical earth is ob- 
tained by multiplying the near field E by the attenuation function, W, as indicated above [1]. 
To provide an illustration of the applicability of the numerical results, the total phase cor- 
rection $( = $0+00 at 100 kc, is shown in figure 10 for the case of a transmitter on land, a 
distance d 2 from the coast. The plotted curves show how <£ varies with d(=d l J [-d 2 ) for various 
values of d 2 . As the receiver moves away from the transmitter, the rate of increase of the 
phase is characteristic of propagation over land. As the coastline is crossed at d=d 2 , the 
phase lag is reduced somewhat and eventually continues to increase at a much slower rate as the 
receiver moves out to sea. This abrupt reduction of the rate of increase of phase lag is known 
as a "recovery effect" and has been verified experimentally [7]. 

2. Concluding Remarks 

The results presented here should be useful in predictions of ground-wave propagation at 
low radiofrequencies over a two-part mixed land-sea path. It should be emphasized that the 
sky waves have not been considered. At ranges greater than 200 km or so, the ionospherically 
reflected waves must be separately accounted for. In any case, the ground wave is omni- 
present even out to 2,000 km and is a substantial contribution to the field. 

It is of interest to note that a formula almost identical to eq (13) of part 1 has been reported 
by Godzinski [8] in Poland at the recent International Radio Consultive Committee meeting in 
Warsaw. His results, however, are not valid for short distances or very low frequencies, where 
the induction and static fields are significant. He has also developed a formula essentially 
equivalent to eq (6) of the present paper, when (d/a) 3 <^<^l. Furutsu [9] has also recently 
investigated the problem of ground-wave propagation over an inhomogeneously conducting 
spherical earth. He has indicated that the basic integral equation can be solved by an iterative 
procedure. The first term of his rather elaborate analysis agrees with eq (7) of the present 
paper. Neither Godzinski or Furutsu present any explicit numerical result for the mixed-path 
attenuation function. 

The results in this paper, as well as in part 1, are not valid for points close to the boundary 
because a first-order stationary-phase evaluation of the integrals has been employed. Further- 
more, the wave reflected from the boundary has been neglected. A somewhat different analyti- 
cal approach has been used [10] to study the field very close to the boundary, which indicates 
the stationary-phase approximation is valid if A and B are at least two wavelengths from the 
boundary or coastline. 
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